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1.  INTRODUCTION 


The  purpose  of  this  report  is  to  explain  a  new  mathematical  method  of  measuring  the  range, 
deflection,  and  height  of  burst  for  a  fired  artillery  shell.  This  new  mathematical  procedure  can  be  used 
as  an  alternative  to  triangulation  (Roberts  1990).  It  can  be  employed  when  there  arc  three  or  more 
tower  readings  of  azimuth  and  elevation  on  a  fntd  projectile. 

In  triangulation,  when  three  or  four  tower  readings  of  azimuth  and  elevation  arc  observed  and  the 
smallest  optimum  triangle  whose  area  is  less  than  50  square  meters  is  found,  range,  deflection,  and 
height  of  burst  are  computed  as  follows.  Range  is  the  distance  from  the  gun  to  center  of  "impact 
point"  in  the  found  triangle.  Deflection  from  line  of  fire  is  computed  as  a  perpendicular  distance  from 
"impact"  to  the  line  of  lire  and  as  an  angle  between  line  of  fire  and  vector  from  weapon  to  "impact 
point."  Altitudes  are  computed  from  each  tower  using  the  associated  elevation  and  tower  height  and 
are  then  averaged  (Roberts  1990). 

The  new  technique  operates  differently.  A  least-squares  method  which  minimizes  angular  error  is 
employed  iteratively  to  find  the  range,  deflection,  and  height  of  burst  point  while  making  use  of  the 
same  data  input  to  triangulation.  The  procedure  is  used  for  all  tower  readings  as  well  as  for  all 
combinations  when  one  tower  reading  is  omitted.  Then,  statistically,  error  computations  and 
associated  burst  points  can  be  compared  to  see  if  there  is  significant  difference.  Thus,  this  method 
provides  a  means  of  detecting  a  bad  tower  azimuth/elevation  observation  in  the  process  of  using  all  the 
information  one  has  access  to. 

This  report  will  first  present  the  algorithms  used  in  this  technique.  Then,  a  discussion  of  the 
algorithms  and  associated  computer  programs  (Appendices  A,  B,  C,  and  D)  will  follow.  The 
conclusion  will  discuss  the  results  of  computer  runs  by  both  triangulation  and  the  new  technique 
applied  to  the  same  observational  data  input. 

2.  ALGORITHMS 

2.1  Least-Squares*  Algorithm  (developed  by  Mr.  Robert  Lieske  and  used  in  subroutine  LSTQR  in 
least-squares  computer  program). 

f  For  more  information  on  least-squares,  see  "Introduction  to  Statistical  Theory,"  by  P.  G.  Hocl,  S.  C.  Pon,  and 
C.  T.  Stone,  1971. 
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X 


(X„  Y„  Z)  -  position  coordinates 
of  burst  point  formed  by 
two  intersecting  tower  azimuth 
readings  and  an  elevation  from 
one  tower. 


Step  1.  For  burst  point  position  (Xjt  Y,,  Z,),  compute  root  mean  square  errors  in  azimuth  (ERMS)  and 
elevation  (ERMSEL)  from  n  tower  reading  locations. 

Let, 

A  AZ*  =  azimuth  angle  (of  segment  from  tower  k  to  burst  point)  -  AZ(I,K)  (azimuth  reading  of  round 
from  tower  k),  and 

A  EL^  =  elevation  angle  (of  segment  from  tower  k  to  burst  point)  -  ALZ(I.K)  (elevation  reading  of 
round  from  tower  k). 

Then, 


ERMS  = 


.1/2 


j,  (AAZt)2 


;  ERMSEL  = 


.1/2 


f  (AElt)2 


(1) 


Step  2.  Find  AX,  AY,  AZ  so  that  the  following  quantities  are  minimized  where  from  a  tower  position 
(X(k),Y(k),Z(k)),  one  has  for  k  =  l,..,n: 

Xk  =  X;  -  X(k); 

Yk=Y,-  Y(k); 

Zk  =  Zi-  Z(k); 


*.i 


dAZ.  dAZ. 

AAZ  -  _1  AY  -  __i 
*  1)AX  UaF 


AY 


dAZt 

3az 


AZ 


(2) 
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N 

*'2  =  E 

Km  1 

Take  the  partial  derivative  of  E2  in  Equation  2  with  respect  to  AX  and  AY  and  take  the  partial 
derivative  of  E'2  in  Equation  3  with  respect  to  AZ.  Set  each  partial  derivative  to  zero  and  rearrange 
terms  to  give  the  following  equations: 


dEL.  dEL.  dEL. 

A EL.  -  AX  -  -  AX  -  _  *  AZ 

*  lAX  dAY  7)AZ 


(3) 


r  \ 

dAZ, 


y  ^AX  j 


■  E 


dAZ, 


v~5ax  ; 


ax+E 


faAzt 


dAZ, 


vw, 


Ay+E 


3az, 


v^AZy 


dAZ, 


^Haxj 


AZ  ; 


(4) 


E^J 


(  \ 

dAZ k 


V3"/ 


X  A 
3AZt 


v3ax, 


/  > 


v9AYy 


kx  +  E 


(  Y 
dAZ, 


K  dAY  j 


AX  +  E 


(  \ 
dAZ, 


ydAZ  , 


(  ^ 

0AZt 


V8AX, 


AZ  ; 
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> 
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w 
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(6) 


Note  the  following  terms  can  be  derived: 


dAZt 

Hax 


X]  +  X2 


dAZk 

Hay 


x2+x2 


3£Lt_  -Xt  Zt 
dA*  (X2+X2+Z2)(X*+X2)1/2 


dAZ. 

_ *  =0 

’  -5az 


(&=l,..,n)  ; 
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MLt_  -Yk  Zk 
^AY  (X2k  +  Y2t+Zl)(X2k+Y2kyn 


6ELk  ~(X2k  +  Y2k)m 

SZz  (xI+yI+zI) 


With  the  above  terms  substituted  into  Equations  4,  5,  6,  one  notices  that  Equations  4  and  5 
become  two  equations  in  the  two  unknowns  AX,  AY,  and  Equation  6  is  one  equation  in  the  unknowns 
AX,  AY,  AZ. 

Solve  Equations  4  and  5  simultaneously  to  obtain  AX,  AY.  Substitute  the  values  AX,  AY  just 
obtained  into  Equation  6  so  the  AZ  can  be  solved  for.  One  has  just  computed  AX,  AY,  AZ  so  that 
Equations  2  and  3  have  been  minimized. 

Step  3.  If  AX,  AY,  AZ  <  .1,  then  the  point  taken  as  the  burst  point  is  (Xt,  Yit  Z;).  Otherwise, 
compute  new  Xit  Yit  Z,  values  as  follows: 

X;  =  X;  -  AX; 

Y,  =  Y,  -  AY; 

Zj  =  Zj  +  AZ; 

and  then  repeat  steps  1,  2,  and  3.  If  the  process  does  not  converge  in  200  iterations,  a  message, 
"DOES  NOT  CONVERGE,"  is  printed  out. 

This  algorithm  can  be  applied  to  the  case  of  all  towers  and  to  the  case  where  each  combination  of 
one  tower  is  omitted. 


2.2  Modified  F-Test  Algorhhrr^With  Outlier  Scheme  -  (used  in  subroutines  MODF  and  OUTLY 
in  least-squares  computer  program). 


*  For  more  information  on  the  modified  F-Test  algorithm,  see  "Introduction  to  Statistical  Analysis,"  by  W.  J. 
Dixon  and  F.  J.  Massey,  1957. 
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Using  algorithm  2.1),  one  obtains  the  following  values: 


ERROR  (J) 
ERRORE(J) 
XCNT(J) 
YCNT(J) 
ZCNT  (J) 


-  Erms  (azimuth) 

-  Erms  (elevation  -  el) 

-  X  coordinate  of  burst  point 

-  Y  coordinate  of  burst  point 

-  Z  coordinate  of  burst  point 


when  omitting  tower  J  (J  =  l,..,n)  and  when  all  towers  are  used  the  following: 

ERROR  (JALL)  -  Erms  (azimuth) 

ERRORE(JALL)  -  Erms  (elevation  -  el) 

XCNT  (JALL)  -  X  coordinate  of  burst  point 
YCNT  (JALL)  -  Y  coordinate  of  burst  point 
ZCNT  (JALL)  -  Z  coordinate  of  burst  point 

Step  1.  Find  smaliest  ERROR(J)  that  is  less  than  ERROR(JALL).  If  any,  conduct  modified  F-tcst 
with  outlier  scheme  as  follows: 


-  ERR0R(K)  . 

*  ~h 


p  =  £  ERROm 0  . 

i u\  n 

K*J 


*♦1 

S  =  Y,  (ERRORS 0  -  x  f  ; 

Km  1 


*♦1 

S'  *  Y  (ERR0R(fC)-7)2  ; 

tf.l 

KtJ 


H«1 


Y 

S 
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Step  2.  Check  to  see  if  Fn+I  derived  from  modified  F-test  is  significant  (i.e.,  is  an  outlier)  using 
n+1  degrees  of  freedom.  Subroutine  OUTLY  in  least-squares  computer  program  has  data  statements 
giving  criterion  for  10%,  5%,  2.5%  and  1%  levels  of  significance  and  up  to  8  degrees  of  freedom 
(Spiegel  1975;  Hoel,  Port,  and  Stone  1971). 

If  ERROR(J)  is  found  to  be  significantly  smaller  then  ERROR(JALL),  then  along  with  range, 
deflection,  burst  height  (HOB),  and  error  values  of  burst  point  found  using  all  towers,  similar  values 
for  burst  position  found  with  tower  J  omitted  are  also  printed  out. 

Step  3.  Find  smallest  ERRORE(J)  that  is  less  than  ERRORE(JALL).  If  any,  conduct  modified 
F-test  with  outlier  scheme  in  a  manner  similar  to  Steps  1  and  2. 

3.  DISCUSSION 

After  running  the  least-squares  computer  program  (Appendix  B)  modified  with  special  write 
statements  on  the  sample  observational  data  input  of  Appendix  C,  the  process  of  convergence  of  the 
least-squares  algorithm  was  investigated.  It  was  found  that  when  using  all  towers  or  when  a  tower 
was  omitted,  convergence  took  place  in  three  or  less  iterations  of  the  subroutine  LSTQR.  Also,  the 
root  mean  square  errors  had  the  tendency  to  reduce  in  the  iteration  process.  This,  of  course,  is 
expected.  When  the  convergence  process  was  begun  at  a  different  initial  burst  point  (Xj.Y;,  Z)  by 
changing  the  combination  of  intersecting  tower  reading  azimuths  and  tower  elevation,  it  ended  at  the 
same  point.  That  is,  although  the  initial  (X,,  Yit  Zj)  of  the  convergent  process  was  changed,  the  final 
(Xj,  Yi,  Z;)  was  invariant  and  produced  the  same  identical  output  in  Appendix  D. 

It  must  be  admitted  that  the  convergent  process  of  the  least-squares  algorithm  was  not  studied  in 
great  detail.  A  statement  of  theorem  dealing  with  convergence  or  divergence  of  the  algorithm  was 
unable  to  be  made  and  proven  by  the  author.  Perhaps  some  pure  mathematician  who  reads  this  report 
will  make  such  an  attempt. 

4.  CONCLUSION 

The  least-squares  computer  program  which  utilizes  the  least-squares  algorithm  is  contained  in 
Appendix  B  with  associated  information  in  Appendix  A.  This  computer  program  along  with  Dennis 
Flaherty’s  triangulation  computer  program  has  been  applied  to  the  same  observational  data  input 
(Appendix  Q,  and  the  results  are  shown  in  Appendices  D  and  E.  Appendices  D  and  E  show  that  the 


6 


least-squares  computer  program  either  with  all  towers  or  with  a  tower  omitted  give  results  in 
deflection,  range,  and  height  of  burst  that  are  similar  to  the  results  of  the  triangulation  program  for 
rounds  1 1  through  17.  Round  18,  however,  is  an  interesting  case.  Depending  on  the  computer 
program,  the  range  values  for  this  round  differ  absolutely  by  some  20  meters.  Which  program  does 
one  believe  in?  The  author  believes  the  results  of  the  least-squares  computer  program  because  it  uses 
all  the  information  and  there  is  no  statistical  reason  for  eliminating  any  of  the  four  tower  azimuth 
readings. 

Note  that  the  least-squares  computer  program  can  incorporate  time  of  flight  observational 
measurements  in  the  same  manner  as  the  triangulation  program  (Roberts  1990).  Also,  the  least-squares 
computer  program  has  been  written  for  the  Aberdeen  Proving  Ground  (APG)  where  azimuths  arc 
measured  clockwise  from  due  south,  and  the  Cartesian  coordinate  system  at  the  gun  has  its  x-axis 
oriented  35°  clockwise  from  due  south,  with  positive  z-axis  going  into  the  ground.  This  program  can 
easily  be  adapted  to  the  coordinate  systems  of  other  proving  grounds.  The  least-squares  computer 
program  is  on  the  VAX  8600  and  VAX  780  computers  in  the  Launch  and  Flight  Division  of  the  U.S. 
Army  Ballistic  Research  Laboratory,  APG,  MD. 
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DEFINITIONS  OF  INPUT  TO  LEAST— SQUARES  COMPUTER  PROGRAM 


XT, YT,ZT 

ITOWER 

IP 

X<7) ,Y(7) ,Z( 
AZDG , AZMG 

-  AZD , AZM 

EZD , EZM 
IDl , ID2 


Position  coordinates  of  T-th  tower  (T  . LE .  6)  in  APG  grid  system 
Tower  number  (1,2,3, 4, 5, 6) 

Any  character  other  than  a  blank 
7)  Position  coordinates  of  weapon,  APG  system 

Line  of  fire  measured  clockwise  from  south,  in  degrees  and 
minutes 

Tower  azimuth  readings  measured  from  south,  in  degrees  and 
minutes 

Tower  elevation  readings  of  burst  height, in  degrees  and  minutes 
Round  identification 
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LISTING  OF  INPUT  TO  LEAST-SQUARES 
All  input  is  punched  into  a  file  using 

I  terns 

XT , YT  #  ZT , I TOWER , I P 

XN,YN,ZN, ITOWER, IP  N  . LE .  6 

Blank 

X(7),Y(7),Z(7) , AZDG , AZMG 

For  N  of  Towers: 

AZD ( I , J ) , AZM , EZD , EZM , IDl , ID2 

If  blank  data  use  -999.  for  both  AZD  and  EZD 


AZD ( K , J ) 
+  999  . 


COMPUTER  PROGRAM 

the  following  format. 
Format 

3Fl0.2,47X,Il, IX ,Al 

3  Fl 0 . 2 , 2 FI  0 . 1 
2F7.2/2F7.2/2x,2A3 
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DESCRIPTION  OF  OUTPUT  TO  LEAST- SQUARES  COMPUTER  PROGRAM 


.  no . 

Deflections 
mils  meters 

Range 

meters 

Hob 

meters 

Erms-Az 

mils 

E  rms-El 
mils 

Towe  r  s 

2A3 

F6 . 1 

F6 . 1 

F8 . 1 

F7.1 

F8 . 3 

F8. 3 

ALL 

2A3 

F6  .1 

F6 . 1 

F8 .1 

F7.1 

F8  .  3* 

F8 . 3 

OMITING 
TOWER  J 

2A3 

F6 . 1 

F6 . 1 

F8 . 1 

F7.1 

F8 . 3 

F8 . 3  * 

OMITTING 
TOWER  J ' 

*  Line 

printed 

only  if 

value  is 

found  signi 

ficant  in 

outlier 

subroutine 
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APPENDIX  B: 

LISTING  OF  LEAST-SQUARES  COMPUTER  PROGRAM 
WRITTEN  FOR  ABERDEEN  PROVING  GROUND,  MD,  BY  NEAL  ROBERTS 
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Intentionally  left  blank. 
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non 


PROGRAM  TO  FIND  RANGE , DEFLECTION , AND  HEIGHT  OF  BURST 
OF  FIRED  ARTILLERY  PROJECTILES  BY  USING  A  LEAST 
SQUARES  TECHNIQUE  WHICH  MINIMIZES  ANGULAR  ERROR. 

IMPLICIT  REAL* 8 ( A-H , O-Z  ) 

DIMENSION  AZ (201,7) ,ALX(201,7) ,AZD(201,7)  ,AZM(201,7) 

DIMENSION  EZD ( 201,7), EZM( 201 , 7),X(7),Y(7),  Z(7) 

DIMENSION  AZMLS ( 201,7),  EZMLS( 201 , 7 ) , I SB ( 7 ) ,CNT( 3 ) , PT( 3 ) 
DIMENSION  RNG ( 7 , 201  )  ,XCNT( 7 ) , YCNT ( 7 )  ,  ZCNT(  7  ) 

DIMENSION  ERROR( 7 ) , ERROREt 7 ) 

DIMENSION  IDl ( 201 ) , ID2 ( 201 ) 

WRITEt  6 , 700  ) 

N=0 

1  READ (5,501)  XT, YT, ZT, ITOWER, IP 
IF(IP.EQ.1H  )  GO  TO  2 
X( ITOWER ) =XT 
Y ( ITOWER )=YT 
Z ( ITOWER ) =ZT 
N  =  N+ 1 
GO  TO  1 

2  READ( 5,502)  X(7),Y(7),Z(7), AZDG , AZMG 
AZ  ( 1 , 7  )  =  (  ( AZDG+AZMG/60. DO )-35. DO >*.017453293 
IRDG=0 
DO  70  J= 1 , N 
K=  1 

DO  70  1=1,200 
GO  TO  (30,60)  K 

30  READ (  5 ,1030)  AZD( I , J ) , AZM( I , J) , EZD ( I , J ) , EZM  (  I ,  J )  ,  IDl ( I ) , ID2 ( I ) 

32  IF( AZD( I , J ) . EQ.-999 .DO )  GO  TO  60 

33  IF( AZD( I , J ) .NE . +999 . DO )  GO  TO  50 

34  IF( IRGD-I )  35,40,40 

35  I RDG= I 
40  K=2 

GO  TO  60 

50  IF(EZD(I,  J)  .EQ.-999.D0)  GO  TO  60 

AZ ( I , J )  =  ( ( AZD ( I , J )+AZM( I , J )/60 .DO ) - 3  5 .DO >*.017453293 
ALX ( I , J ) = ( EZD ( I , J ) +EZM( I, J)/60. DO)*. 017453293 
GO  TO  70 

60  AZ(I,J)=99999.D0 
ALX ( I , J ) =-999 . DO 
70  CONTINUE 

IRDG=IRDG-1 

80  DO  81  1=1, I RDG 

81  AZ (1,7) =AZ (1,7)  ' 

C  WRITE  HEADINGS 

WRITEt  6 , 701 ) 

WRITEf  6 , 702 ) 

DO  280  1=1,1 RDG 
NHR=1 
J  =  1 

DO  90  K= 1  ,  N 

IF(AZ(I,K) .EQ. 99999. DO)  GO  TO  90 

NHR=NHR+1 

ISB ( J ) =K 

J-J  +  l 

90  CONTINUE 
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GO  TO  (274,110,120,130,130,130,130) ,NHR 
274  WRITE(6, 1000)  IDl ( I ) , ID2 ( I  ) 

GO  TO  280 
110  I  SB ( 2 ) =7 

120  CALL  POINT(AZ , I , ISB,X, Y , PT) 

CNT ( 1 )=PT( 1 ) 

CNT ( 2 ) =PT ( 2 ) 

CALL  RANGE (X(7),Y(7),AZ(1,7), CNT , RN , DEFL , DEFLM ) 

CNT ( 1 ) =PT ( 1 ) 

CNT ( 2 ) =PT ( 2 ) 

I F ( I  SB ( 2 ) . NE . 7 )  GO  TO  121 
KK=I SB ( 1  ) 

CALL  RANGE ( X ( KK  )  , Y ( KK ) , AZ ( 1 , 7 ) , CNT , RNG ( KK , I ), DEFL , DEFLM ) 

DZ=RNG ( KK , I ) *TAN ( ALX ( I , KK ) ) 

ZN=Z ( KK ) +DZ 
GO  TO  122 

121  KKl  =  I  SB ( 1  ) 

KK 2= I  SB ( 2 ) 

CALL  RANGE(X(KK1),Y(KK1)  ,AZ( 1,7), CNT, RNG(KKl, I), DEFL, DEFLM) 

CALL  RANGE ( X ( KK2  )  ,Y(KK2)  ,  AZ  ( 1 , 7 ) , CNT , RNG ( KK2 , I ) , DEFL , DE FLM ) 
DZl=RNG{ KKl , I ) *TAN ( ALX ( I , KKl ) ) 

DZ2=RNG ( KK2 , I ) *TAN ( ALX ( I , KK2 ) ) 

ZNl=Z(KKl)+DZl 
ZN2=Z( KK2 )+DZ2 
ZN=(ZNl+ZN2)/2. 

122  WRITE( 6 , 1002 )  IDl(I),ID2(I), DEFLM, DEFL,  RN ,  ZN 
GO  TO  280 

130  NT=NHR-1 

CALL  POINT ( AZ , I , ISB,X,Y,PT) 

ISB(NHR)=ISB(NT)+1 
CNT( 2 )=PT( 2 ) 

CNT( 1 )=PT(1 ) 

KK= I SB ( 1  ) 

CALL  RANGE(X( KK) , Y ( KK ) , AZ ( 1 , 7 ) , CNT, RNG (  KK,  I )  ,  DEFL,  DEFLM) 

DZ=RNG( KK, I ) *TAN(ALX( I , KK) ) 

PT(3)=-(Z(KK)+DZ) 

DO  140  IT=1 , NHR 
CNT ( 3 ) =PT ( 3 ) 

CNT ( 2 ) =PT( 2 ) 

CNT ( 1 ) =PT ( 1 ) 

I TER= I  SB ( IT) 

CALL  LSTQR(CNT, ERMS, ITER, I , I SB , AZ , NT, X , Y , Z , ERMSEL , ALX , IQ , IDl , ID2 ) 
IF( IQ. EQ. 1 )  GO  TO  280 
XCNT( ITER ) =CNT ( 1  ) 

YCNT ( I TER ) =CNT ( 2  ) 

ZCNT( ITER ) =CNT ( 3 ) 

ERROR ( ITER ) =ERMS 
ERRORE( ITER)=ERMSEL 
140  CONTINUE 
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JALL=ISB(NHR) 

ERMS- ERROR ( JALL ) 

CNT ( 1 )=XCNT( JALL ) 

CNT(2)=YCNT( JALL) 

CALL  RANGE  (  X  (  7  )  ,  Y  (  7  )  ,  AZ  (1,7)  ,  CNT ,  RN,  DEFL ,  DEFLM  ) 

ERMSE=ERRORE( JALL) 

HOB=-ZCNT( JALL) 

WRITE (  6 ,1003)  I  Dl  (  I ) , I D2 ( I ) , DEFLM , DEFL , RN , HOB , ERMS , ERMSE 
IF(ERROR( JALL)  . EQ . 0 . .AND. ERRORE( JALL)  .  EQ  .  0  .  )  GO  TO  280 
IF(ERROR( JALL)  . EQ . 0  )  GO  TO  150 

CALL  MODF(  ERROR,  I  SB  ,  NHR ,  NT ,  I  TOWER ,  VARl ,  VAR ,  1 1 1  ) 

IF(III.EQ.O)  GO  TO  150 

CALL  OUTLYf VARl , VAR, SIG, ITEST,NHR) 

IF( ITEST. EQ. 0 )  GO  TO  150 

J=ITOWER 

CNT ( 1 ) =XCNT ( J  ) 

CNT ( 2 ) =YCNT { J ) 

ERMS=ERROR ( J ) 

ERMSE=ERRORE ( J ) 

HOB  =  -ZCNT( J  ) 

CALL  RANGE  (X(  7)  ,  Y (  7  )  ,  AZ  (  1 , 7  )  ,  CNT ,  RN,  DEFL  ,  DEFLM  ) 

WRITE ( 6 ,1004)  IDl ( I ) , ID2 ( I ) , DEFLM ,  DEFL  ,  RN,  HOB , ERMS , ERMSE, J  ,  SIG 
150  IF(ERRORE( JALL) . EQ . 0)  GO  TO  280 

CALL  MODF(  ERRORE  ,  I  SB,  NHR,  NT,  I  TOWER,  VARl ,  VAR,  III  ) 

IF( III . EQ. 0  )  GO  TO  280 

CALL  OUTLY ( VARl , VAR, SIG, ITEST, NHR) 

IF( ITEST. EQ. 0)  GO  TO  280 

J=ITOWER 

CNT ( 1 ) =XCNT ( J ) 

CNT ( 2 ) =YCNT ( J ) 

ERMS=ERROR( J ) 

ERMSE  =  ERRORE ( J  ) 

HOB=-ZCNT( J ) 

CALL  RANGE  (  X  (  7  )  ,  Y  (  7  )  ,  AZ  {  1 , 7  )  ,  CNT ,  RN ,  DEFL  ,  DEFLM  ) 

WRITE ( 6 ,1005)  IDl ( I ) , ID2( I ) , DEFLM ,  DEFL ,  RN,  HOB , ERMS , ERMSE,  J,  SIG 
280  CONTINUE 
STOP 

501  FORMAT ( 3F10.2, 47X, Il,lX,Al) 

502  FORMAT ( 3F10.2,2F10.1) 

700  FORMAT ( ' 1 ' ) 

701  FORMAT (  'O'  , ' RD . NO .  '  ,T12, 'DEFLECTIONS'  ,T30,  'RANGE'  ,T40,  'HOB'  ,T55, 'E 
1RMS-AZ '  , T6  5 ,  'ERMS-EL'  ,T80, 'TOWERS' ) 

702  FORMATf 1H  , Til ,' MILS ', T19 , 'METERS' ,T39, 'METERS  '  ,  T57 , ' MILS ' , T67 , ' MI 
1LS  *  ) 

1000  FORMAT! '0' ,A3,A3) 

1002  FORMAT!  '0'  ,  2A3 , T9 , F6 . 1 ,  T19  ,  F6 . 1 ,  T28 , F8 . 1 , T38 , F7 . 1 ) 

1003  FORMAT!  '0'  ,  2A3 , T9 , F6 . 1 , T19 , F6 . 1 , T28 , F8 . 1 , T38  ,  F7 . 1 , T54 , F8 . 3 , T6 4 , F3 . 
13 ,T82 , 'ALL'  ) 

1004  FORMAT! 1H  , 2A3 , T9 , F6 . 1 ,  T19  ,  F6 . 1 ,  T28 , F8 . 1 , T38  ,  F7 . 1 , T54 , F8 . 3 , T6 2 , '*' 
1,T64,F8.3,T80, ' OMITING  TOWER  ',11,'  (SIG.  AT  ',F4.2,'  LEVEL)') 

1005  FORMAT! 1H  ,  2A3 , T9 , F6 . 1 , T19 , F6 . 1 , T28 , F8 . 1 , T38  ,  F7 . 1 ,  T54 , F8 . 3 , T6 4 , F8 . 
13, T72, '*' ,T80, 'OMITING  TOWER  ',11,'  (SIG.  AT  ',F4.2,'  LEVEL)') 

1030  FORMAT! 2F7 . 2 , 2F7 . 2 , 2X, 2A3 ) 

END 


21 


SUBROUTINE  POI NT ( AZ , I , I  SB , X , Y , PT ) 

IMPLICIT  REAL *8(A-H,0-Z) 

DIMENSION  AZ(201,7),ISB(7),X(7),Y(7),PT(3) 

J I SC= I  SB ( 1  ) 

KI SC= I  SB ( 2  ) 

A=TAN(AZ( I, JISC)  ) 

B  =  TAN ( AZ ( I , RISC)  ) 

PT ( 1 )  =  ( X ( JI SC ) *A-Y ( JI SC ) +Y ( RISC ) -X ( KI SC ) *B  )  /  (  A-B  ) 
PT(  2 )  =  ( PT( 1 )-X( RISC)  ) *B  +  Y( RISC) 

RETURN 

END 
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SUBROUTINE  RANGE ( X, Y, AZ , CNT , RNG , DEFL  ,  DEFLM  ) 
IMPLICIT  REAL  *  8 ( A-H , 0-2 ) 

DIMENSION  CNT ( 3  ) 

RNG=CNT ( 1 )-X 
DEFL=CNT ( 2 ) -Y 

CALL  ATNSR(DEFL, RNG, BETA, 1) 

DEFLM=BETA-AZ 

10  I F ( CEFLM+ i .14x59265}  30,30,20 
20  I F ( DEFLM-3 . 14159265  )  50,40,40 
30  DEFLM-DEFLM+6 . 28318531 
GO  TO  10 

40  DEFLM=DEFLM-6 . 28318531 
GO  TO  10 

50  RNG=RNG/COS( BETA) 

DEFL=RNG* S IN ( DEFLM) 

DEFLM=1018 . 59164*DEFLM 

RETURN 

END 
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SUBROUTINE  ATNSR ( A , B , C , L ) 
IMPLICIT  REAL *8 ( A-H , O-Z ) 
P-3. 14159265 
C-ATAN ( A/B ) 

I F ( A )  10,40,40 
10  I F ( B )  20,30,30 

20  GO  TO  ( 21 , 22  )  ,  L 

21  C=P+C 

GO  TO  50 

22  C=C-2 . * P 
GO  TO  50 

30  GO  TO  ( 31 , 50)  ,L 

31  C=2 . *  P+C 
GO  TO  50 

40  I F( B )  20,50,50 
50  RETURN 
END  . 


SUBROUTINE  LSTQR ( CNT , ERMS , I TER , I , I  SB , AZ , NT ,  X  ,  Y , Z , ERMS EL , ALX , I Q , I Dl 
1 , ID2 ) 

IMPLICIT  REAL *8 ( A-H , O-Z  ) 

DIMENSION  AZ ( 201 , 7  )  ,ALX( 201 , 7  )  ,CNT( 3 ),ISB(7),X(7),Y(7),Z(7) 
DIMENSION  DAZ ( 6 ) ,DAZX(6) ,DAZY(6) ,RNG(7,201) 

DIMENSION  DEL ( 6 ) ,DELX( 6  )  ,  DELY ( 6 ) , DELZ ( 6 ) ,IDl(201),ID2(201) 

I COUNT=l 
IQ=0 

183  I COUNT= I  COUNT + 1 

I F ( I COUNT . NE . 2  0 1 )  GO  TO  10 
WRITE( 6,1)  IDl ( I ) , ID2 ( I ) 

1  FORMAT( IX, 2A3 , 5X, ' LSTQR  SUBROUTINE  DID  NOT  CONVERGE  IN  200  ITERATI 
IONS' ) 

IQ=1 
GO  TO  9 

10  DO  190  11=1, NT 
J=ISB( II ) 

IF( J.EQ.ITER)  GO  TO  190 
YC=CNT( 2 ) -Y ( J ) 

XC=C NT(1)-X(J) 

CALL  ATNSR(YC, XC, ANGLE, 1) 

DAZ (II) =ANGLE-AZ ( I , J ) 

184  IF(DAZ( II )+3 .14159265)  186,186,185 

185  IF(DAZ( II )-3. 14159265)  189,188,188 

186  DAZ ( II )=DAZ( II )+6 . 28318531 
GO  TO  184 

188  DAZ( II )=DAZ( II )-6 . 28318531 
GO  TO  184 

189  CONTINUE 

DAZX ( II )=-YC/(XC**2+YC**2 ) 

DAZY ( II )=XC/(XC**2+YC**2 ) 

190  CONTINUE 
N=0 

SERMS  =  0  . 

DO  191  11=1, NT 
J=ISB( II  ) 

IF( J.EQ. ITER)  GO  TO  191 
N=N+1 

SERMS  =  S ERMS  +  DAZ ( 1 1  )  *  *  2 

191  CONTINUE 

ERMS  =  1018 . 5916  4  *SQRT ( SERMS/FLOAT( N ) ) 

TM1  =  0  . 

TM2  =  0  . 

TM3=0 . 

TM4=0 . 

TM5=0 . 

TM6  =  0  . 

DO  192  11=1, NT 
J= ISB (II) 

I F ( J.EQ. ITER)  GO  TO  192 
TM1=TM1+DAZ( II ) *DAZX( II ) 

TM2=TM2+DAZX( 1 1 ) *  *  2 
TM3=TM3+DAZX( II ) *DAZY( II ) 

TM4-TM4+DAZ (II) *DAZY( II ) 

TM5  =  TM5  +  DAZX( II ) *DAZY(  1 1  ) 

TM6=TM6+DAZY ( 1 1 ) *  *  2 

192  CONTINUE 
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DY=( TM1*TM5-TM4*TM2 )/( TM3 *TM5-TM6 *TM2 ) 

DX  = ( TM1-TM3  *DY )/TM2 
DO  200  11=1, NT 
J  =  I  SB ( II ) 

IF( J . EQ. ITER)  GO  TO  200 
YC=CNT (  2  )  -  Y  (  J  ) 

XC=CNT ( 1 ) -X ( J ) 

ZC=CNT( 3  )  -rZ  (  J  ) 

CALL  RANGE ( X ( J ) , Y ( J  )  ,  AZ  ( 1 , 7  )  ,  CNT , RNG ( J , I ) , DEFL , DEFLM ) 
ANGLE=ATAN ( ZC/RNG ( J , I ) ) 

DEL ( II ) =-ANGLE-ALX( I , J ) 

194  IF(DEL(II)+3. 14159265)  196,196,195 

195  I F ( DEL (II)-3. 14159265)  199,198,198 

196  DEL( II )=DEL( II )+6 . 28318531 
GO  TO  194 

198  DEL( II )=DEL( II )-6 . 28318531 
GO  TO  194 

199  CONTINUE 

Dl= ( XC*  *2+YC**2  +  ZC**2 ) 

D2=SQRT( XC**2+YC**2 ) 

DELX ( II ) = ( -XC*ZC )/( D1 *D2 ) 

DELY (II)=(-YC*ZC)/(Dl*D2) 

DELZ ( II )=D2/D1 

200  CONTINUE 
N=0 

SERMS  =  0  . 

DO  201  11=1, NT 
J=ISB (II) 

I F ( J . EQ. ITER)  GO  TO  201 
N=N+1 

SERMS=SERMS+DEL ( II ) **2 

201  CONTINUE 

ERMSEL=1018 . 59 1 6 4 *SQRT ( SERMS/FLOAT  (  N )  ) 

TM1  =  0 . 

TM2  =  0  . 

TM3  =  0  . 

TM4  =  0  . 

DO  202  11=1, NT 
J=ISB ( II  ) 

I F ( J . EQ. ITER)  GO  TO  202 
TM1  =  TM1  +  DEL ( II ) *DELZ( II  ) 

TM2=TM2+DELX( II ) *DELZ ( II ) 

TM3=TM3  +  DELY( II  ) *DELZ ( II ) 

TM4=TM4+DELZ ( 1 1 ) *  *  2 

202  CONTINUE 

DZ»(TM1-TM2*DX-TM3*DY)/TM4 

I F ( ABS ( DX )  . LE .  .1 .AND . ABS ( DY ) . LE .  .  1  .AND . ABS ( DZ ) . LE  .  .  1 )  GO  TO  9 
CNT ( 2 ) =CNT ( 2 ) — DY 
CNT ( 1 )=CNT( 1 )-DX 
CNT (3)=CNT(3)+DZ 
GO  TO  183 
9  RETURN 
END 


26 


SUBROUTINE  OUTLY ( ERMS , ERROR , S IG , I TEST ,  NT  ) 

IMPLICIT  REAL *8 ( A-H , O-Z ) 

DIMENSION  F99 ( 8 ) , F975( 8 ) , F95( 8 ) , F90 ( 8 ) 

DATA  F99/.0,  .0,.  0001,. 0100,. 0442,.  0928,.  .1447,-194  8/ 
DATA  F97  5/.0  , .0  ,  .0007  ,  .02  48,  .0808, .  14  53, .2066, .2616/ 
DATA  F9S/.  0  , .  0  ,  .0027  ,  .0494  ,  .1270  , .2032, .2696,  .3261/ 
DATA  F90/.0  , .0, .0109,  .097  5, .1984  , .2826, .3503  ,  .4050/ 
F=ERMS/ERROR 

IF(F.LE.F99(NT) )  GO  TO  10 
IF( F . LE . F97 5 ( NT ) )  GO  TO  20 
I F ( F . LE . F9 5 ( NT ) )  GO  TO  30 
I F ( F . LE . F9 0 ( NT )  )  GO  TO  40 
ITEST=0 
RETURN 
10  SIG=  .  99 
I TEST= 1 
RETURN 
20  SIG= . 97 5 
ITEST=1 
RETURN 
30  SIG= . 95 
ITEST= 1 
RETURN 
40  SIG= . 90 
I TEST= 1 
RETURN 
END 
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SUBROUTINE  MODF( ERROR, I  SB , NHR , NT , I TOWER , VARl , VAR , III ) 
IMPLICIT  REAL* 8 { A-H , O-Z ) 

DIMENSION  ERROR ( 7  )  ,  I  SB ( 7 ) 

TEMP= ERROR ( NHR ) 

Xl=ERROR ( NHR ) 

DO  160  I T= 1 , NT 
J=ISB( IT) 

ERMS=ERROR ( J ) 

I F ( ERMS . GE . TEMP )  GO  TO  160 
ITOWER= J 
Xl  =  ERROR ( J  ) 

TEMP=ERROR ( J ) 

160  CONTINUE 

I F ( Xl . GE . ERROR ( NHR ) )  GO  TO  169 
SUM=0 . DO 
DO  162  I T= 1 , NHR 
J=ISB( IT) 

SUM=SUM+ERROR ( J ) 

162  CONTINUE 

AVG= SUM/FLOAT ( NHR ) 

SUM=0 .DO 
DO  164  IT=1 , NHR 
J  =  I  SB ( IT) 

I F ( J.EQ. ITOWER)  GO  TO  164 
SUM  =  SUM+ERROR ( J  ) 

164  CONTINUE 

AVG1 =SUM/FLOAT ( NT ) 

VAR=0 . DO 
DO  166  I T= 1 , NHR 
J=ISB( IT) 

VAR=VAR+ ( ERROR ( J ) -AVG ) *  *  2 
166  CONTINUE 
VAR1=0 . DO 
DO  168  I T= 1 , NHR 
J=ISB ( IT) 

IF (J.EQ. ITOWER)  GO  TO  168 
VAR1=VAR1+ ( ERROR( J ) -AVG1 ) **2 

168  CONTINUE 
III-l 

GO  TO  170 

169  111=0 

170  RETURN 
END 
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APPENDIX  C: 

SAMPLE  INPUT  TO  LEAST-SQUARES  COMPUTER  PROGRAM 
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Intentionally  left  blank. 
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3657  , 

.607 

12801 .626 

32.13 

3868  . 

.530 

13870.487 

13  .  3 

3754  , 

.335 

13094 . 383 

4.77 

2883  , 

.79 

11998  .  31 

12.13 

3812  , 

.  273 

12475.803 

41.0 

34  . 

22  . 

18  . 

51  . 

11 

34  . 

23  . 

17  . 

45  . 

12 

34. 

25. 

18. 

46  . 

13 

34  . 

18  . 

17  . 

41  . 

14 

34  . 

11  . 

5. 

35. 

15 

34  . 

17  . 

5  . 

09  . 

16 

34  . 

05. 

5. 

35. 

17 

34  . 

+  999  . 

17. 

5. 

08  . 

18 

11  . 

08. 

19. 

07. 

11 

9. 

54  . 

18. 

01 . 

12 

11  . 

13  . 

19  . 

05. 

13 

S  . 

42  . 

17  . 

42. 

14 

10  . 

52  . 

6  . 

09. 

15 

11 . 

19  . 

6  . 

37  . 

16 

10. 

36  . 

6. 

01. 

17 

11 . 

+  999  . 

03  . 

5. 

32. 

18 

26. 

57. 

20  . 

04  . 

11 

27. 

28  . 

18  . 

50. 

12 

27  . 

52. 

19  . 

45. 

13 

27. 

24  . 

18  . 

40  . 

14 

27  . 

42  . 

6. 

45. 

15 

27  . 

57. 

5. 

55. 

16 

27  . 

32. 

6  . 

35. 

17 

27  . 

+  999. 

45  . 

5. 

50  . 

18 

47. 

37. 

14  . 

48. 

11 

47  . 

55. 

13. 

38  . 

12 

47  . 

35. 

14  . 

56  . 

13 

48  . 

03  . 

13  . 

41  . 

14 

47  . 

26. 

4  . 

35. 

15 

47  . 

24  . 

4  . 

19. 

16 

47. 

33  . 

4  . 

39  . 

17 

47  . 

+  999. 

35. 

4  . 

17  . 

18 

31 


Intentionally  left  blank. 
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APPENDIX  D: 

SAMPLE  OUTPUT  FROM  LEAST-SQUARES  COMPUTER  PROGRAM 


33 


Intentionally  left  blank. 
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RD.NO.  DEFLECTIONS  RANGE  HOB  ERMS-AZ  ERMS-EL  TOWERS 

MILS  METERS  METERS  MILS  MILS 


r3  *-3 

►3  J 

ca  ca 

ca 

ca 

ca 

ca  ca 

Ea  & 

Ea 

Ea 

Ea 

Ea  Ea 

►a  j 

>j 

►3  J 

in  as 

o 

in 

in 

in  <j\ 

©  Ch 

on 

cr> 

cn 

on  o> 

o  o 

o 

o 

o 

o  o 

Eh  Eh 

E-< 

E* 

Eh 

Eh  Eh 

<  < 

< 

< 

< 

<  < 

d  d 

d 

d 

d 

d  d 

HH  M 

M 

h-t 

w 

HH  >— t 

to  to 

to 

to 
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APPENDIX  E: 

SUMMARIZED  DATA  OUPUT  TO  TRIANGULATION  COMPUTER  PROGRAM  USING  SAMPLE 

DATA  INPUT  OF  APPENDIX  C 
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RD.NO. 

DEFLECTIONS 

RANGE 

MILS 

METERS 

METERS 

11 

4 . 3 

10.9 

2551.5 

12 

7.7 

18.3 

2433.0 

13 

4 . 4 

11.1 

2562.7 

14 

9.4 

22.3 

2407.5 

15 

1.0 

2.6 

2543.3 

16 

0.8 

1.9 

2587.1 

17 

0.2 

0.5 

2520.9 

18 

3.5 

8.6 

2530.4 

AVERAGE  HOB 
METERS 

952.9 

853.9 

952.9 
840.2 
302.1 
292.8 
297.5 
273.0 
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No  Of  No  Of 

Copies  Organization  Copies  Organization 


2  Administrator 

Defense  Technical  Info  Center 
ATTN:  DTIC-DDA 
Cameron  Station 
Alexandria,  VA  22304-6145 

1  HQDA  (SARD-TR) 

WASH  DC  20310-0001 

1  Commander 

US  Army  Materiel  Command 
ATTN:  AMCDRA-ST 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333-0001 


1  Commander 

US  Army  Missile  Command 
ATTN:  AMSMI-RD-CS-R  (DOC) 

Redstone  Arsenal,  AL  35898-5010 

1  Commander 

US  Army  Tank- Automotive  Command 
ATTN:  AMSTA-TSL  (Technical  Library) 
Warren,  Ml  48397-5000 

1  Director 

US  Army  TRADOC  Analysis  Command 
ATTN:  ATRC-WSR 

White  Sands  Missile  Range,  NM  88002-5502 


1  Commander 

US  Army  Laboratory  Command 
ATTN:  AMSLC-DL 
Adelphi,  MD  20783-1145 


(cim».  only)  i  Commandant 

US  Army  Infantry  School 
ATTN:  ATSH-CD  (Security  Mgr.) 
Fort  Benning,  GA  31905-5660 


2  Commander  <uncias*.  only)  i  Commandant 


US  Army,  ARDEC 

ATTN:  SMCAR-IMI-I 

Picatinny  Ars  ?nal,  NJ  07806-5000 

2  Commander 

US  Army,  ARDEC 

ATTN:  SMCAR-TDC 

Picatinny  Arsenal,  NJ  07806-5000 

1  Director 

Benet  Weapons  Laboratory 
US  Army,  ARDEC 
ATTN:  SMCAR-CCB-TL 
Watervliet,  NY  12189-4050 

1  Commander 

US  Army  Armament,  Munitions 
and  Chemical  Command 
ATTN:  SMCAR-ESP-L 
Rock  Island,  IL  61299-5000 

1  Director 

US  Army  Aviation  Research 
and  Technology  Activity 
ATTN:  SAVRT-R  (Library) 

M/S  219-3 

Ames  Research  Center 
Moffett  Field,  CA  94035-1000 


US  Army  Infantry  School 
ATTN:  ATSH-CD-CSO-OR 
Fort  Benning,  GA  31905-5660 

1  Air  Force  Armament  Laboratory 
ATTN:  AFATUDLODL 

Eglin  AFB,  FL  32542-5000 

Aberdeen  Proving  Ground 

2  Dir,  USAMSAA 

ATTN:  AMXSY-D 

AMXSY-MP,  H.  Cohen 
1  Cdr,  USATECOM 

ATTN:  AMSTE-TD 

3  Cdr,  CRDEC,  AMCCOM 

ATTN:  SMCCR-RSP-A 
SMCCR-MU 
SMCCR-MSI 

1  Dir,  VLAMO 

ATTN:  AMSLC-VL-D 
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No.  of 

Conies  Organization 

1  Director 

HQ,  TRAC  RPD 

ATTN:  ATRC-RP 

Fort  Monroe,  VA  23651-5143 

1  Director 
TRAC-FLVN 
ATTN:  ATRC 

Fort  Leavenworth,  KS  66027-5200 

2  Commander 

US  Anny,  ARDEC 
ATTN:  SMCAR-CAWS-S, 

Mr.  R.  DcKlcinc 
Mr.  D.  Griggs 

Picalinny  Arsenal,  NJ  07806-5000 

1  OPM  Nuclear 

ATTN:  AMCPM-NUC 
Picalinny  Arsenal,  NJ  07806-5000 

7  Commander 

US  Army,  ARDEC 
ATTN:  SMCAR-AET, 

Mr.  F.  Sccrbo 
Mr.  J.  Bera 

ATTN:  SMCAR-AET-A, 

Mr.  R.  Kline 
Mr.  H.  Hudgins 
ATTN:  SMCAR-FSA, 

Mr.  F.  Brody 
Mr.  R.  Kantenwein 
ATTN:  SMCAR-FSS, 

Mr.  J.  Brooks 

Picalinny  Arsenal,  NJ  07806-5000 


No.  of 

Conics  Organization 

1  Commander 

US  Army  Dugway  Proving  Ground 
ATTN:  STEDP-MT 

Mr.  G.  C.  Travers 
Dugway,  UT  84022 

1  Commander 

US  Army  Yuma  Proving  Ground 
ATTN:  STEYP-MTW 
Yuma,  AZ  85365-9103 

1  Headquarters 

US  Marine  Corps 
ATTN:  Code  LMW/30 
Washington,  DC  20380 

1  Director 

Sandia  National  Laboratories 
ATTN:  Mr.  A.  Hodapp 
Division  1631 
Albuquerque,  NM  87185 

1  Commander 

Naval  Surface  Warfare  Center, 
Aerodynamics  Branch, 

K-24,  Building  402-12 
ATTN:  Dr.  W.  Yanta 
White  Oak  Laboratory 
Silver  Spring,  MD  20910 

1  Director 

Lawrence  Livermore  National  Laboratory 
ATTN:  Mail  Code  L-35 
Mr.  T.  Morgan 
PO  Box  808 
Livermore,  CA  94550 


2  Commandant 

US  Army  Field  Artillery  School 
ATTN:  ATSF-CCM 
ATTN:  ATSF-GD 
Fort  Sill,  OK  73503 

1  Director 

US  Army  Field  Artillery  Board 
ATTN:  ATZR-BDW 
Fort  Sill,  OK  73503 
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No.  of 

Conics  Organization 


Aberdeen  Proving  Ground 

2  Director,  USAMSAA 
ATTN:  AMXSY-RA, 

Mr.  R.  Scungio 
AMXSY-GS, 

Mr.  B.  King 

1  Commander,  USATECOM 
ATTN:  AMSTE-TE-F, 

Mr.  W.  Vomocil 

l  PM-SMOKE,  Bldg.  324 
ATTN:  AMCPM-SMK-M, 
Mr.  J.  Callahan 

1  Director,  USAHEL 
ATTN:  SLCHE-FT 

2  Commander,  USACSTA 
ATTN:  STECS-AS-H 

STECS-EN-B 


Intentionally  left  blank. 


USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes. 
Your  comments/answers  to  the  items/questions  below  will  aid  us  in  our  efforts. 

1.  BRL  Report  Number _ BRL -MR-3894 _ _  Date  of  Report  February  1991 

2.  Date  Report  Received  _ _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest 

for  which  the  report  will  be  used.)  _ _ _ _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source 
of  ideas,  etc.) _ _ _ _ _ 


if 

5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hour's  or  dollars 
saved,  operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate 
changes  to  organization,  technical  content,  format,  etc.) _ _ _ 


- * - 

- '  J* 

'•  *  _ 

Name 


CURRENT  Organization 

ADDRESS  _ 

Address 


City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  New  or  Co.rci'f 
Address  in  Block  6  above  and  the  Old  or  Incorrect  address  below. 


Name 


OLD  Organization 

ADDRESS  _ 

Address 


City,  State,  Zip  Code 


(Remove  this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  mail.) 


FOLD  HERE 


Department  of  the  Army 

Director 

t'.S.  Army  Ballistic  Research  Laboratory 

ATTN:  SLCBR-DDT 

Aberdeen  Proving  Ground,  MD  210i  -5066 

OFFICIAL  BUSINESS 


BUSINESS  REPLY  MAIL 

FIRST  CLASS  PERMIT  No  0001,  APG,  MD 


POSTAGE  WILL  BE  PAID  BY  ADDRESSEE 


Director 

U.S.  Army  Ballistic  Research  Laboratory 

ATTN:  SLCBR-DD-T 

Aberdeen  Proving  Ground,  MD  21005-9989 


FOLD  HERE 


